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Abstract 

We formulate the infrared regularization of Becher and Leutwyler in a form analogous to our 
recently proposed extended on-mass-shell renormalization. In our formulation, IR regularization 
can be applied straightforwardly to multi-loop diagrams with an arbitrary number of particles 
with arbitrary masses. 
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I. INTRODUCTION 



Starting with Weinberg's fundamental work on phenomenological Lagrangians ^ 
became possible to systematically calculate corrections to the soft-pion results obtained 
within the framework of current algebra [2]. The corresponding effective field theory 
(EFT) — chiral perturbation theory (ChPT) — has been very successful in describing the 
strong interactions at low energies (for a recent review see, e.g., Ref. Q). In the mesonic 
sector, the combination of standard dimensional regularization (DR) and the modified 
minimal subtraction scheme of ChPT (MS) led to a straightforward correspondence be- 
tween the loop expansion and the chiral expansion in terms of momenta and quark masses 
at a fixed ratio 0, 11 • The one-baryon sector proved to be more complicated || . In par- 
ticular, using the same combination of DR and MS as in mesonic ChPT, higher-order 
loops contribute in lower chiral orders and therefore the correspondence between the loop 
expansion and the chiral expansion seems to be lost (see Fig. 2 of Ref. One solution 
to this power-counting problem was given in the framework of heavy-baryon chiral per- 
turbation theory (HBChPT) |7(, and most of the recent calculations have been performed 
within this approach jHEJ. While successful in many cases, HBChPT destroys the ana- 
lytic structure in part of the low-energy region. Several methods have been suggested to 
reconcile power counting with the constraints of analyticity in a manifestly relativistic 
approach 11 , 12, [r| [13, Ell- The one most widely used is the infrared (IR) regu- 
larization of Ref. ll| by Becher and Leutwyler. A possible generalization to multi-loop 
diagrams has been suggested in Ref. 

In the present paper we provide a formulation of the IR regularization of Becher and 
Leutwyler in a form analogous to the extended on-mass-shell (EOMS) renormalization of 
Ref. [15]. As a result of the reformulation, IR regularization is straightforwardly applicable 
to multi-loop diagrams as well as to diagrams involving several fermion lines and/or 
resonances. 



II. COMPARISON OF IR REGULARIZATION AND EOMS RENORMALIZA- 
TION 

In order to reformulate the infrared regularization of Becher and Leutwyler [ll| in a 
form analogous to the EOMS renormalization of Ref. [15j , we consider as an example the 
characteristic, dimensionally regularized, one-loop integral of the fermion self-energy, 

_ . r d n k 1 

nA-P, p _ p y _ m 2 + [fc2 _ M 2 + 10 + ] ' O 

where n denotes the number of space-time dimensions. The masses m and M refer to 
the nucleon mass in the chiral limit and the lowest-order pion mass, respectively. Using 
the standard power counting of Refs. |l7], EH we assign the order Q n ~ 3 to the integral 
ijVTr- Here, Q collectively denotes small expansion parameters such as the pion mass or 
small external momenta. (Note that I^n satisfies power counting only after subtraction 
10). 



To implement the IR regularization and to compare with the EOMS renormalization 
scheme we use the Feynman parametrization formula 

L=[ l ^ (2 ) 

ab Jo [az + bil-z)} 2 ' V ' 
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with a = (k—p) 2 — m 2 + iO + and b = k 2 — M 2 +iO + , interchange the order of integrations, 
perform the integration over loop momenta k, and obtain 



1 Ntt\ 



-P, 0) 



1 



r(2-n/2) / dz [A{z)\ 



,(n/2)-2 



(3) 



where 



= -p 2 (l - 2)2 + m 2 z + M 2 (l 



2) - iO" 



In the approach of Becher and Leutwyler, the integral In-k is divided into the IR singular 
part / and the remainder R, = I + R, defined as 



R 



(4tt)™/ 2 
1 



r(2-n/2) f°° dz[A(z)] (n/2) - 2 



dz [A(z)\ 



(n/2)-2 



(4) 



(5) 



( 4vr )n/2 

In this decomposition, for noninteger n the integral I is proportional to a noninteger 
power of the pion mass (~ M n ~ 3 ) and thus satisfies the power counting. On the other 
hand, the remainder R does not satisfy the power counting but, for arbitrary n, contains 
non-negative powers of small parameters and is thus absorbed into an infinite number of 
counterterms. Divergent parts of I are also absorbed in an infinite number of counter- 
terms. 

In our EOMS renormalization scheme 2| we apply a conventional renormalization 
prescription which allows us to identify the terms which we subtract from a given integral 
without calculating the integral beforehand. In essence we work with a modified inte- 
grand which is obtained from the original integrand by subtracting a suitable number of 
counterterms. To find the subtraction terms we consider the series 



{p 2 -m 2 ) l {M 2 y 




1 



[(k - p) 2 -m 2 + i0+] [k 2 — M 2 + i0 + ] 

1 



p 2 =m 2 ,M 2 =0 



(k 2 - 2k-p + i0+) (k 2 + i0~ 



+ (P 2 



m 



2m 2 (k 2 
1 



2k ■ p + i0- 



(k 2 -2k-p + i0+y (k 2 + i0+)_ 



(k 2 - 2k-p + i0+) (k 2 + i0+y 
1 1 
2m 2 (k 2 - 2k ■ p + i6+) (k 2 + i0~ 



+ 



(6) 



where [. . .] p 2 =m 2 means that we consider the coefficients of (p 2 — m 2 ) l (M 2 y only for four- 
momenta p^ which satisfy the on-mass-shell condition. Although the coefficients still 
depend on the direction of after integration of this series with respect to the loop 



momenta k and evaluation of the resulting coefficients for p 



m 



the integrated series 



is a function of p 2 only. We subtract from Eq. those terms of the expansion of Eq. (jUJ) 
which violate the power counting. These terms are analytic in the small parameters and 
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do not contain infrared divergences. For the given example we only need to subtract the 
first term of the expansion of Eq. © • 

We note that integrating Eq. © term by term reproduces the expansion of R of Eq. (jSJ) 



in M 2 and p 2 



m 



This can be checked by explicitly integrating the first few coefficients 



of the expansion of Eq. (EJ 
expansion of R of Ref. |llj 



we indeed see that they coincide with the coefficients of the 



R 



m n " 4 r(2 - n/2) 
(47r)"/ 2 (n - 3) 



m 



in 



6)(P 2 



[n 



3)M 2 



2m 2 



Am 4 (n - 5) 



2m 2 (n — 5) 



+ 



(7) 



A more straightforward and transparent way of obtaining the IR regular part R is to 
rewrite I^-n using the Feynman (or Schwinger) parameterization, integrate over loop mo- 
menta, expand the resulting integrand (of the integration over parameters) in a Taylor se- 
ries of Lorentz-invariant small expansion parameters (small masses and Lorentz-invariant 
combinations of external momenta and large masses), and, finally, interchange summation 
and integration: / dx — > S J dx. As is shown in the next section, the above observation 
is correct in general, i.e., by expanding the integrand of any integral with an arbitrary 
number of nucleon and pion denominators in small parameters and interchanging summa- 
tion and integration, one reproduces the expansion of the IR regular part of the integral. 1 



III. GENERAL CASE 

Let us consider the general one-loop scalar integral corresponding to diagrams with 
one fermion line and an arbitrary number of pion and fermion propagators: 



I N ...ir...{Pl, ...,?!,...)=« J 

where 



d n k 1 



(2tt)» b 1 ---b l a 1 



bj = (k + pj) 2 - m 2 + i0 + , 
Oi = (k + qif - M 2 + i0 + . 

Tensor integrals are reduced to the scalar integrals of Eq. (jHJ) in the standard fashion ^| . 

Following Ref. 0] we apply the infrared regularization to the integral of Eq. (jHJ). We 
start by combining all meson propagators using the formula 

l ( d \ {m - 1] r\ r 1 X 

dxx • • • / dx m -i-r . (9) 



a\---a m \dM 2 ) Jo Jo A 

The numerator X is given by 



X 



1 for m = 2, 

X2(x 3 ) 2 ■ ■ ■ (x m _i) m_2 for m > 2, 



1 Note the important difference with Ref. [lo| . where the expansion of the integrand with subsequent 
interchange of integration and summation reproduces the chiral expansion of the power-counting pre- 
serving part. 
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and the denominator A is given by the recursive expression 

A = A m , 
Ai = at, 

A p+1 = x p A p + (1 - x p )a p+1 (p = 1, . . . , m - 1). 
The result for A is of the form 

A= (k + q) 2 - A + iO + , (10) 

where the constant term A is of order Q 2 , and q is a linear combination of external 
momenta and is of order Q l . 

Analogously we combine the nucleon propagators 

1 ( d r 1 f 1 Y 

dyi-" / dyi-x—. (11) 



bi • • • bi \ dm 2 J Jo Jo B 

The numerator Y is given by 



Y 



1 for I = 2, 

y 2 (y 3 ) 2 ■ ■ ■ (yi-i) 1 - 2 for / > 2, 



and the denominator B is given by the recursive expression 

B = B u 
Bi = h, 

Bp+i = y P B p + (1 - y p )b p+1 (p = 1, . . . , I - 1). 
The result for I? reads 

B = (k + P) 2 — B + i0 + , (12) 

where P is a linear combination of external momenta, P 2 = m 2 + 0(Q) and B = m 2 + 
0{Q). 

Next we combine the denominators A and 5 using 

1 r 1 dz 



AB Jo [{I - z)A + zB] 2 
and obtain for the integral of Eq. (JHJ) 

*(^) (^) Io dz Io dVl '"Io dyi - 1 

rl rl r \ 

x / dxi--- l dx m -\YX I - — r t>. (13) 

io Jo J (2tt)" [(1 - + z£] 2 V ; 

Substituting A and £> from Eqs. (JTUJ) and ()12jl in Eq. (JT3J), evaluating the derivatives, and 
shifting k — » k — Pz — q(l — z), we obtain 

i(I + m - 1)1 [ it - ,)-» | ^ • • • /; ^ O- / ^L [k2 J {z)]l+m , (14) 



where 



f(z) = P 2 z 2 -[P 2 -B)z + A(l -z)- (q 2 -2P-q) z{\ - z) - iO + . 
Finally, the integration of Eq. (|14|) over k yields 

u \n/2 ni + m-n/2) dzz l -\\-z) m - x \ dy x - ■ • / dx m ^YX[f(z)] (n ^- l - m . 
(An) n ' z Jo Jo Jo 

(15) 

To apply the IR regularization we rewrite the z integration as 

1 roo roo 

dz ■ ■ ■ = dz • • • — I dz ■ ■ ■ . 

o Jo h 

The result of the first integration is identified as the IR singular part and of the second 
as the IR regular part. In the IR regular part one can expand the integrand in small 
momenta and masses and interchange summation and integration [ill ]. This leads to 
integrals over z of the type 

/•oo 

Ii= / dzz n+i , (16) 



where i is an integer number. These 7, are multiplied by (further) integrals over Xj 
and yk which do not depend on n. The integrals of Eq. are calculated by analytic 
continuation from the domain of n in which they converge, i.e. 



z n+i+l 



n + i + 1 



n + i + 1 



' (17) 



On the other hand, if we expand the integrand in Eq. ()15j) in small momenta and masses 
and interchange summation and integration, we obtain exactly the same expansion as for 
the IR regular part of the IR regularization with the only difference that instead of the 
integrals Jj of Eq. (J 16)) we now have 2 



./, = - I dzz n+i . 







Calculating these integrals by analytical continuation from the domain of n in which they 
converge, we obtain: 

. - i 

1 



Ji = - 



n + i + 1 



n + i + 1 



(19) 



Comparing Eqs. (|T7|) and ()19j) we see that the expansion of the integrand in Eq. ()15j) with 
subsequent interchange of summation and integration exactly reproduces the result of the 
IR regular part of the loop integral. Next we observe that, if we expand the integrand of 
Eq. ()14J) in small parameters and interchange summation and integration over k, we ob- 
tain exactly the same result as by expanding the integrand in Eq. (|To^ in small parameters 
with subsequent interchange of summation and integration over Feynman parameters. We 
further note that the result of the expansion of the integrand of Eq. (|14|) in small parame- 
ters with subsequent interchange of summation and integration over k coincides with the 



The minus sign relative to Eq. I|ltj[l stems from the definition of R as — dz ■ 
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series which is obtained when we formally expand the integrand of the original integral 
in small parameters, using a formula analogous to Eq. (jUJ), interchange summation and 
integration, and rewrite the integrals of the obtained series in Feynman parametrization. 
We thus conclude that the IR regular part of the original integral can be obtained by 
expanding the integrand in small parameters and interchanging summation and integra- 
tion over loop momenta. In practical calculations of the IR regular parts of loop integrals 
it is convenient to reduce the loop integrals to integrals over (Feynman/Schwinger) pa- 
rameters, expand the integrand in Lorentz-invariant small expansion parameters (small 
masses and Lorentz-invariant combinations of external momenta and large masses), and 
interchange integration and summation. 



IV. APPLICATIONS 

As a check and application of our formulation of the IR regularization we have explicitly 
verified for all integrals of pion-nucleon scattering of Ref . |20j (to the order which is needed 
for the accuracy of calculations of that work) that, by expanding the integrands in small 
parameters and changing the order of summation and loop integration, one reproduces 
the IR regular parts of these integrals. Although the original integrals of BChPT do 
not contain infrared divergences, the IR regular parts as well as the IR singular parts 
separately contain such divergences. 3 In the approach of Becher and Leutwyler these 
divergences of both parts are absorbed in counterterms. (In fact they exactly cancel each 
other and hence do not give any contributions in couterterms.) In our formulation the 
IR regularized integrals are obtained by subtracting the IR regular parts, from which 
the IR divergences are removed beforehand, from the full expressions of the integrals. 
Clearly our expressions of the IR regularized integrals coincide with the results of the 
Becher-Leutwyler approach. 

It is straightforward to apply our formulation of IR regularization to diagrams with 
multiple nucleon lines. We have checked that our approach reproduces the results of 
Ref. [21] for diagrams with two nucleon propagators. As an illustration let us consider 
the following integral: 4 

d n k 1 



/d k 
(2?T [(* + Pi) 2 - m? + i0+] [{k - P 2 ) 2 -m 2 + i0+] [k 2 — M 2 + z0+] 

(20) 

Using Fe ynm an parametrization and performing loop momenta integration one can write 
I NNw as El 

WPi, -P* 0) = r(3 ~ n/2) f Q dz z £ dw[C(w, z) - z0+] W2M , (21) 



where 

C(w,z) = {l-z)M 2 +zm 2 -z 2 {l-w 2 ) 



(P 1 + P 2 ) 2 z{l-z){P 2 + P 2 ) wz{l - z){P 2 - P 2 



3 Note that, using dimensional regularization, IR divergences are also parametrized as l/(n — 4) poles. 

4 Our notations differ from those of Ref. 21]. 
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Following Ref. [2l| we define the IR regular part of the integral Inntt(Pi, —P2, 0) as 



Rnntt(Pi, —P2, 0) 



1 r(3-n/2) 



(4 vr )n/2 



dz z 



dw + 



dw 



C(w,z)-i0' 



(n/2)-3 



(22) 



To calculate Rnntt we expand the integrand of Eq. (|2^j) in powers of M , 4m — (Pi+P 2 ) , 
P 2 — m 2 , and P 2 — m 2 and interchange integration and summation 2l] . Doing so we obtain 
a series, the coefficients of which are proportional to the integrals 



dz ( z 2)W)-* z 1+i 



dw + 



dw 



(w 2 )^- 3 w^ 



where i and j are integers. Again, the integrals Jy are calculated by analytical continua- 
tion from the domain of n in which they converge, leading to 



-IV 



(n -A + i)(n - 5 + j)' 



(23) 



On the other hand, in our approach we identify the IR regular part of Inn-k by expanding 
the integrand in Eq. (|2*Tj) in powers of small parameters [M 2 , Am 2 — (P± + P2) 2 , P\ — rn 2 , 
and P 2 



m 



and interchanging summation and integration over Feynman parameters. 
This leads to exactly the same expansion that we obtained above for R NN7T (Pi, —P 2 ,0), 
but instead of the integrals 1^ we now have 



J, 



dz(z 



2^(n/2)-3 z l+i 



dw{w 



which we calculate by analytically continuing from the domain of n in which they converge 



J, 



I.I 



(n - 4 + i)(n - 5 + j)' 



(24) 



Clearly, in analogy to the one-nucleon sector our formulation of IR regularization repro- 
duces the results of Ref. 2l( for diagrams involving two nucleon lines. 

Recently, we have shown that, within the EOMS renormalization scheme, one 
can set up a consistent power counting in the effective field theory with (axial) vector 
mesons included explicitly. Analogously we could apply the IR regularization in our 
formulation. When treating vector mesons in the antisymmetric tensor field representation 
and analyzing the diagrams contributing to the electromagnetic form factors of the nucleon 
up to and including C(g 4 ), we observe that in Ref. all relevant loop diagrams have 
actually been taken into account. This is due to the fact that the integrals involving only 
vector meson and nucleon propagators vanish in IR regularization. 

Finally, we have also applied our formulation of the IR regularization to the integral 
considered in Ref. 24]| in the context of including the A resonance: 



'OttI 



-p, 0) 



d n k 



(27r) n [(k - p) 2 + i0+] [k 2 — M 2 + z0+] ' 



For p 2 ^> M 2 the chiral dimension of this integral should be Q 
calculation of the integral 1^ results in 

M n-4 r(2 - n/2)r(n/2 - 1) 



n-2 



'CM 



-p,0) 



(4tt) 



n/2 



r(n/2) 



1,2 



n/2; n/2; 



p 2 + i0 H 
M 2 



(25) 

An explicit 
') , (26) 
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where F(a, b; c; z) is the hypergeometric function [25|. For p 2 > M 2 we rewrite Eq. (126)) 
as 

M n-2 r[l-n/2l / . . M 2 \ 

!oA-P, 0) = fA W9 1 ' ^ F [1,2 — n/2; n/2- ] 



( 47r )n/2 p 2 y ' / ' / ' p 2 + y 

(_ p 2 _ i0 +)(n/2)-2 r(2-n/2)r(ra/2-l)r(7i/2-l) / M 2 \"~ 3 

(4vr)«/ 2 r(n - 2) \ ~ p 2 + i0+ J 

(27) 

Analogously to Ref. 0] we identify the first term in Eq. (|27p . which, for noninteger values 
of n, is proportional to a noninteger power of M, as the IR singular part and the second 
term as the IR regular part. The IR singular part satisfies the power counting and would 
generate an infinite number of terms if the function multiplying M n ~ 2 were expanded in 
powers of M 2 . This differs from the result of Ref. 24l . where only the first term of such 



an expansion was identified as the IR singular part of I 0n (—p, 0). 

In analogy to the self-energy integral considered above, it is straightforward to check 
explicitly that, if one expands the integrand of Eq. (|25|) in powers of M 2 and interchanges 
integration and summation, one exactly reproduces the expansion of the second (regular) 
term of Eq. (|27|) in powers of M 2 /p 2 . 

The integral Ishl~ P, 0) has an imaginary part for p 2 > M 2 in both definitions, ours 
and that of Ref. [24] [provided that in this work the same boundary condition has tacitly 
been assumed as in Eq. 1)25)1]. which is included in the IR regular part. This imaginary 
part is given by 

' [l-M 2 /p 2 ) 



16tt 

in n = 4 dimensions and violates the power counting. Therefore, although the regular 
part is analytic in M 2 and consequently its real part can be absorbed by counterterms 
of the Lagrangian, the imaginary part cannot be altered. As a result there exists no 
subtraction scheme within which the renormalized version of Ion{—p, 0) would satisfy the 
power counting. However, from this observation one should not draw the conclusion that 
there is no consistent power counting in a manifestly Lorentz-invariant formulation of 
BChPT with spin 3/2 particles included explicitly. Rather, as already pointed out in Ref. 
23|, the integral l 0n (—p,0) occurs when the spin 3/2 particle propagator is decomposed 



using projection operators and the apparent puzzle disappears once the results for this 
decomposition are put together. 



V. SUMMARY AND CONCLUSIONS 

We have reformulated the IR regularization of Becher and Leutwyler in a form 
analogous to our EOMS renormalization scheme of Ref. (l5T |. Within this (new) formu- 
lation the subtraction terms are found by expanding the integrands of loop integrals in 
powers of small parameters (small masses and Lorentz-invariant combinations of external 
momenta and large masses) and subsequently exchanging the order of integration and 
summation. Isolating the infrared divergences from these terms and subtracting them 
from the original relativistic loop integral one obtains the IR regularized expression of 
the integral. Within our formulation of IR regularization one does not necessarily need 
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to use dimensional regularization. The renormalized results are regularization-scheme in- 
dependent, provided that the infrared divergences are properly handled, and one could 
thus also use, e.g., Pauli-Villars regularization. One advantage of the new formulation of 
IR regularization is that it can be easily applied to diagrams with an arbitrary number of 
propagators with various masses (e.g., resonances) and/or diagrams with several fermion 
lines as well as to multi-loop diagrams. 
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